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We apply the background field method to calculate the moments of the pion two-particles 
twist-3 distribution amplitude (DA) </> p (£) in QCD sum rules. In this paper, we do not use the 

O ! 

equation of motion for the quarks inside the pion since they are not on shell and introduce a 

new parameter to be determined. We get the parameter TOq w 1.30GeV in this approach. 
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obtain its approximate expression which can be determined by its first few moments. 
PACS numbers: 13.20.He 11.55.Hx 



INTRODUCTION 



O 

The perturbative QCD theory has been applied extensively to inclusive and exclusive processes 



during the past decades. The structure functions in inclusive processes, Fi(x, Q ), can be solved by 



matching the theoretical predictions with the experimental values at different Q 2 regions. Unlike 

<D _ 
^ . the structure functions, the distributive amplitude in exclusive processes, <pi(x,Q ), can not easy 



to be measured directly from experiments. The distribution amplitudes in hadrons are the key 
ingredients and provide the universal non-perturbative inputs for many exclusive processes. Al- 
though distribution amplitudes satisfy the renormalization group equation or the QCD evolution 
equation at short distances [l], their solutions depend on the initial conditions (j>i(x, Qq), which are 



determined by the non-perturbative theory. 

n 

In Ref. [2], QCD sum rules were used to study the leading-twist distribution amplitude of the 
pion at the first time. They worked out the first few moments of the twist-2 distribution amplitude 
of the pion and used them to restrict its shape. After that, much work have been done to study 
the leading twist distribution amplitude of pion and their conclusion is that the experimental data 
favors the asymptotic (j>as(Q = f (1 ~ £ 2 ) which is the solution of the QCD evolution equation as 
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Q 2 -»■ oo . 

In the past twenty years, the studies on leading-twist distribution amplitude of the pion also 
had been done by many papers in Ref. DDDD . From the naive point of view, the contribution 
of higher twist distribution amplitudes is suppressed by a factor 1 /Q 2 in large momentum transfer 
processes. In many cases, we can consider solely the contribution from leading-twist distribution 
amplitude, but not involve the impact of higher-twist distribution amplitudes. However the leading 
twist contribution to explain the data for the exclusive processes at the present energy region 
(Q 2 > 4GeV 2 ) leave a lot of space to discuss the contributions from higher twist terms. For 
example, the contribution of twist-3 distribution amplitudes to the pion electromagnetic form 
factor, although it is suppressed by the factor 1/Q 2 , is comparable with and even larger than the 
contribution from leading-twist distribution amplitude of the pion at intermediate energy region 
of Q 2 being 2-40 GeV 2 in Ref Q|. In the semileptonic B meson decays processes, the B — > tt 
transition form factor were calculated in light cone QCD sum rules in the literature The 
B — > tt transition form factor in Ref [(| contain twist-2, twist-3 and twist-4 contributions. From the 
numerical results of Ref. [y] we can see that the contribution of the twist-3 distribution amplitudes 
is about 30-50% and the contribution of the twist-4 distribution amplitudes is about 5% to B — > tt 
transition form factor. Although one can avoid to calculate the twist-3 contribution by choosing an 
adequate chiral current correlator [2j to make the contribution of the twist-3 distribution amplitudes 
in the B — * tt transition form factor vanish automatically, one can't ignore their contributions once 
and for all. On the other hand, we can get the sum rules for the B —* tt transition form factor 
fg^iq 2 ) which is determined by the twist-3 distribution amplitudes solely if choosing an adequate 
chiral correlator 8]. It means that the contribution from the twist-3 distribution amplitudes to the 
fg^iq 2 ) has the same order of magnitude as that from the leading twist distribution amplitude. 
Therefore one has to study the twist-3 distribution amplitudes <j> p (t;), 0o-(£)> (/WC^i) an d higher 
twist distribution amplitudes of the pion and other meson. 

The pion twist-3 distribution amplitudes have been studied in Ref. i], Q] in the chiral limit, 
based on the techniques of nonlocal product expansion and conformal expansion and including 



corrections in the meson-mass. Ref. 



10] studied the vector meson twist-3 distribution amplitudes 



in details. However they employ the equations of motion of on-shell quarks in the meson to get 
two relations among two-particles twist-3 distribution amplitudes of the pion, <^>p(£) and an d 
three-particle twist-3 distribution amplitude cj)^ (oti) of the pion. Then, they took the moments 
of three-particles twist-3 distribution amplitude calculated in QCD sum rules and obtained the 
approximate forms for three twist-3 distribution amplitudes of the pion. The question is whether 
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one can use the equation of motion of the quark inside the meson since the quarks are not on-shell. 
The relations between 4>p(£),<ficr(£) and ^^(ai) can not be satisfied exactly in this point of view. 

In this paper, we do not apply the quark equation of motion and calculate the moments of the 
twist-3 distribution amplitudes <f) p (£) directly in QCD sum rules approach. We employ the QCD 



in the background field 







to calculate the relevant correlator. The (ft p (£) is defined as follow: 



(0|d» 75 exp {ig J* dx^A u(-z)\n+(q)) = m p f^ j d^)e^ + ■■■, (1) 

where the dots refer to those higher twist terms omitted here, and tuq is a parameter which is 
to be determined in our calculation. This parameter has nothing to do with the masses of the 
quarks and/or meson. If we assume the expansion of <fr p (£) in the series in Gegenbauer polynomials 

1/2 

C n (£) one can obtain the approximate expressions of the twist-3 distribution amplitudes of the 
pion which can be determined by their moments. 

For another twist-3 distribution amplitude <^ CT (£) > some words are in order. We found that, 
results calculated directly by QCD sum rules method depend on the Borel parameter seriously. 
That is to say, there is no way to make the 7r-contribution dominate. So we do not consider this 
case in this paper. 



121], axial currents in a correlator 



Finallyit needs to be pointed out that, just as shown in Ref. 
would be coupled to objects like instantons. And this may cause some complication in the calcula- 
tion and make the results unreliable. We will not explore their influences in this paper, some more 
work will be done in later paper. 

This paper is organized as follows: The calculation of two-point correlation function is presented 
in Sec. II. We calculate moments of 4> p {u) in Sec. III. Sec. IV is devoted to build a possible model 
on the twist-3 distribution amplitude. Sec. V is reserved for discussion and summary. 

II. CALCULATION OF TWO-POINT CORRELATION FUNCTION 

In order to calculate two-point correlation function, we apply the background field approach 
which were expressed in Ref. explicitly. They called the QCD in background fields. In that 
framework, the non-vanishing vacuum condensates are described by the classical fields, while the 
corresponding quantum fields are quantized in the Furry representation and the physical states are 
built upon the physical QCD vacuum through the action of creation operators. One can derive the 
propagators of quarks and gluons as a perturbative series in gauge invariant form. The propagators 
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are the solution of equations, 

(i^(d^-igA^-m)S F (x,0) = <5 4 (x), 
( gix „(D p DO) AB + 2f ABC G%)S™(x,0) = 8 AD giMT 5\x). 

To the aim of calculation of first few moments in this paper, we need only consider following terms 
in the two propagators respectively, 

where a,b = 1,2,3, = -if ABC G^(A, B, C = 1,2,..., 8), G% = d^A A - d v A A ^ 

9f ABC A B A^ is the classical field strength, A A is the classical field of the gluon. 
We define a two-point correlation function, 

n ri r 2 (g) = i J d 4 xe^ x (0 | T{ jTl (x)jr 2 (0)} I 0), (4) 

where T± and T\ are Dirac matrices, and the vacuum state is the physical one. 

For the light quark current, the OPE of two-point correlation function is as follow: 



i J d A xe^(0\T{j ri (x) 3r2 (0)}\0) 



= 4%(q 2 )I + 4%(q 2 )(0 | | 0) + C^VW | ^G 2 | 0) + • • • , (5) 

where we ignore the higher power of light quark mass m q and only reserve the linear m q term. 
Cr 1 r 2 ' s i n Eq.(6) are Wilson coefficients which can be worked out in perturbative QCD theory. 
For calculating the moments of <^> p (£) , we take as usual the two-point correlation functions: 

i J d'xe^iO | r{jf n) (x,z)jf f (0)} | 0) = (z ■ q) 2n I 2nfi (q 2 ), (6) 

where the currents are defined as 

ji 2n) (x,z)=d(x) l5 (iz-*D) 2n u(x), (7) 

with D = dfj, — i2gA c ^ J T a . The relation between the light quark currents j^ n ^ (x, z) and the twist-3 
distribute amplitude of the pion will be given in the next section. 



On the one hand, one can calcu 



Following the familiar procedure 



ate the correlation functions in QCD framework perturbatively. 



3] 



by applying operator expansion in the background field ap- 



proach, corresponding Feynman diagrams in Fig.l, with condensates up to dimension-6 and the 
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perturbative contribution part to the lowest order, we have 

T 2n,0 ( _ n 2, _ 3 1 „2,„ f -Q 2 \ _ 1n-l(m u + m d ){4)i)) 



+ — A q 2 Q4— gj-(21 + 8n(n + l)). (8) 



On the other hand, one can express the two-point correlation functions in hadronic represen- 
tation by inserting the complete intermediate states with the same quantum numbers as those of 
corresponding currents which is j 2n in this case. By isolating the pole term of the lowest pseu- 
doscalar pion, the hadronic representation can be obtained in the following, 



:Jd^e^(0\T{i n) (x,z)jP + (0)}\0) 



(0 I jt\z) I ir)fr 1 jf t(Q) | 0) | y (0 | jf n \z) | ir^Trf | j<°> + (0) | 0) 
q 2 -ml 2-, q 2 _ m 2 ff 

where n H s are higher exciting hadronic states. 



III. SUM RULES FOR THE MOMENTS OF TWIST-3 DISTRIBUTION AMPLITUDE 

Twist-3 distribution amplitude <fi p (£) of the pion is given as the leading terms by gauge invariant 
elements of nonlocal hadronic matrices in Eqs.(l). Furthermore, the elements of nonlocal hadronic 
matrices can be expanded near light cone z 2 = 0, 



, ,. z N OO 

(0|J(z) 75 exp I ig / dx^A^ \ u(- z)\tt + (q)) = ]T 
I J-z J „_ n 



n! 

n=0 



(—) n - <— ► 

(0|J(0) 7 5(^ D) n U (0)\ir + (q)). (10) 



Substituting Eq.(14) into Eq.(l), we have the following relations, 

(0\jt l) (z)\n)=f n m p (e p n )(z-qr (11) 
where the moments of </> p (£) is defined as, 

& = l f/ n Mm- (12) 

Therefore the spectrum can be taken as approximately, 

lml 2n '°(q) = v5{q 2 - m 2 )K) 2 /. 2 <C> ■ <0 + vr A ^-^{q 2 - % ), (13) 

where s n is the threshold parameter and should be set to a value between m 2 and its first exciting 
state it'. 
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To get the sum rules for moments, we employ dispersion relation for I(Q 2 ), 



and assume the contribution of the lowest state ir dominant in the sum rule, making a Borel 
transformation upon above dispersion relation, we can get an improved one 



l — 2 Jds e-^Im/f %) = L M I 2n '°(-Q 2 ), 



ttM 2 

where M is the parameter used in the Borel transformation with M = Q 2 /n fixed, 



(14) 



L 



M 



lim 



AQT( 



d 



y 



QV^oo (n - 1)! x ^ ' y dQ 2 ' 

Inserting Eq.(13) and Eq.(8) into the left and right hand side of Eq.(14) respectively, we obtain 
the sum rule for moments of </> p (£), 

3 1 s _jlzl\ 2n - 1 (m u + nid) (V> tp) 



\Sp > Up) (mlfff 



8vr 2 2n + 1 1 ' ~~ ( ' + 
2n + 3(^G 2 ) 16vr 



le m 2 



24 A/ 4 + *•(»+!)) 



2 M 4 
M 6 



(15) 



As we can see, the first term in Eq.(15) is the pure perturbative results with condensate con- 
tributions omitted completely. 



IV. A POSSIBLE MODEL OF THE PION TWIST-3 DISTRIBUTION AMPLITUDE 

Before proceeding further we need to fix input parameters. For the condensate parameters, we 
take as usual, 

{m u + m d ){^) ~ -1.7 x W^GeV 4 , 
(^L G i) ~ 0.014GW 4 , 

7T 

(^alV^) 2 ^ 1-83 x 10" 4 Gey 6 . (16) 

The parameter s n in Eq.(15) should be chosen to make moments and tuq most stable against M 2 
in a certain range. This can be done by considering the ratios of (£,p^ n+1 ^) to (£~ n ) . The advantage 
is that the ratios depend only on s n and M 2 , but not TOq and This make us determine s„- 

and M 2 on a more general ground at first. The results are shown in Fig. 2. It is shown from these 
figures that s w affect the results a little and M 2 should be taken, at least, 1.2GeV 2 . The key point 
to avoid the concept of on-shell equations of motion is that the introduced tUq appearing in the 
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sum rule of moments. And they can be obtained through requiring normalization of the zeroth 
moments((£p) = 1). With these points kept in mind, we can show the results as (M 2 = 1.5— 2GeV 2 ) 

s v = 1.7GeV 2 ml = 1.24 - IMGeV (£ 2 ) = 0.340 - 0.356 {0 = 0.167 - 0.210, 

S7T = l.6GeV 2 ml = 1.19 - l.WGeV = 0.340 - 0.359 (0 = 0.164 - 0.211, (17) 

Sir = l.SGeV 2 ml = 1.14 - 1.2AGeV ($) = 0.341 - 0.361 = 0.160 - 0.212. 

The above numerical results show that the dependence of ml for Borel parameter is about 10%, 
the dependence of (£, 2 ) for Borel parameter is about 5% and the dependence of for Borel 
parameter is about 25% as M 2 is in the region we chose. It is obvious that the fourth moment of 
is more unreliable. The above numerical results also show that the parameter mjj is smaller 



Q3 



when the equation of motion of 



than m 2 /(m u + m^) ~ 1.78 — 3.92GeV which was used in 
on-shell quarks are employed. 

In principle, infinite number of moments are needed to determine twist-3 distribution amplitude 
of the pion completely, that is beyond one's capability by now in this method. Also we have 

neglected contributions from higher dimension condensates in the QCD calculation and higher 
exciting resonances in the hadronic representation. Thus the prediction for more moments is not 
reliable. We will take the above three moments into consideration in the following analysis and 
expand twist-3 distribution amplitudes p (£)of the pion in Gegenbauer polynomials to first few 
terms, 

<f> p (0 =A p + B p C l 2 '\0 + C p C 4 1/2 (0, (18) 

1/2 

where C n (£) are Gegenbauer polynomials. A p ,B p ,C p are three coefficients to be determined 
through the moments given in previous section. We obtain the two group of equations, 



~ £ dtt 2n (A p + B p C l 2 /2 (0 + C p C 4 1/2 (0) = (£f ) n = 0, 1, 2. 



(19) 



Solving the above equations, it is easy to find these coefficients in twist-3 distribution amplitude 

A p = (0 , B p = -§«£> - 3<3» , C, = f(3<#>-30<#> + 36<#>). (20) 

Substituting the value Eq.(17) into Eq.(20), the <j> p (t;) can be obtained. (when (s n = 1.7GeV 2 ), we 
have, A p ~ 1,B„ ~ (0.05) - (0.17), C p ~ (-1.52) - (-0.371). ) We can compare Eq.(18) with the 
results of Ref . 9j , 

MO = 1 + 30i? C 2 1/2 (0 + ^(4u, 2)0 - - 2u lfi ) cl /2 (0, (21) 
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where R = niu r ^-J nd y 21 and the parameters w's can be found from the moments of (j)^ in QCD 
sum rules. Comparison Eq.(21) with our results is expressed in Fig. 3. The distribution amplitude 
depicted in the dashed line contains contributions from the second and the fourth moment. But 
it should be pointed out that the fourth moment calculated in this paper is less reliable than the 
second one. The dotted line depicts the one only contains the contribution from the second moment 
of 4>p(0- Finally, Our result is closed to the asymptotic form, 4>p S (£) = 1- 



V. SUMMARY AND DISCUSSION 



In applying QCD to exclusive processes, form factor and decay amplitude are determined by 
the convolution of a perturbative hard scattering amplitude with a non-perturbative distribution 
amplitude of given twist. The study of the leading twist and non-leading twist distribution ampli- 
tudes is a crucial task for predicting exclusive processes more precisely. The asymptotic behavior 
of exclusive processes is determined by the leading twist distribution amplitude only as Q 2 is very 
large. One has to consider the contributions from the non-leading twist at the present energy 
regions. Many people have claimed that the contributions from the twist-3 distribution amplitudes 
of the pion make a sizeable corrections to the exclusive processes. Therefore it is helpful to study 
the twist-3 distribution amplitudes of the pion in different approaches in order to understand their 
behaviors. 

In this paper, we apply the background field approach to calculate the OPE of the two-point 
correlation functions related to the moments of twist-3 distribution amplitudes of the pion. With 
the aid of the improved dispersion relations, we obtain the sum rules of the moments of (j> p (Q- 
Inputting all sorts of parameters required in sum rules with the contribution of condensate to 
dimension-6, we worked out first three moments of 4> p (£,) ■ It is shown that the main contribution 
to moments is from pertubative term, the contributions from condensate terms are relatively small, 
In the previous approach 0, Q], it was based on the techniques of nonlocal operator product 
expansion and conformal expansion. They applied the equations of motion to impose relations 
between the higher twist distribution amplitudes and 4>3 n . We emphasize that the equations of 
motion can not be satisfied since the quarks are not on mass-shell within hadrons. 

To avoid imposing equations of motion of on-shell quark, it is not felicitous, we introduce a 
new parameter mj in the definitions of moments of (ft p (£) ■ Because of the emergence of the new 
parameter, we can not determined all moments if we don't know . The strategy we take is 
to require normalization condition = 1, and then we obtain rrig and other moments. In this 
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paper, we have made some approximations, such as the lowest pole dominates and the higher 
dimension condensates are negligible, also we haven't calculated the radiative corrections. Based 
on these approximations, we obtain the parameter m.Q ~ 1.30GeV ±0.06GeV which is smaller than 
m^/{m u + nid)- However, it is closed to the phenomenological value Q|. 

Using the values of the parameter ttiq and the moments, the twist-3 distribution amplitude 
4> p {0 expanded in Gegenbauer polynomial can be obtained , 

MO = 1 + 0-137 C 2 1/2 (£) - 0.721 cl /2 (0- 
Its behavior is closed to the shapes of (fip S (0 in Ref.jj]]. 
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FIG. 1: Several types of Feynman diagrams calculated in this paper. The gluon condensations are included in 
the Fermion propagators implicitly and the quark condensations are depicted as the crosses in the diagrams. 
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M 2 GeV 2 

FIG. 2: The ratios of (£,p ) to (^") versus Borel parameter M. The upper curves are the ratios with n=l 
and the lower curves are the ratios with n=0, while the solid, the dashed and the dotted curves correspond 
to three threshold values: (a) = 1.7GeV 2 ; (b) s w = 1.6GeV 2 (c) = 1.5GeV 2 , respectively. 
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FIG. 3: The twist-3 distribution amplitude </> p (£). The solid-curve is the one of Ref.|a ,the dashed curve 
corresponds to Eq.(18) and the dotted curve corresponds to the one only contained the second moment in 
this paper. 



